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Abstract. A multichannel analysis of cosmic ray electron and positron spectra and of the 
diffuse synchrotron emission of the Galaxy is performed by using the DRAGON code. This 
study is aimed at probing the interstellar electron source spectrum down to E < 1 GeV and 
at constraining several propagation parameters. We find that above 4 GeV the e~ source 
spectrum is compatible with a power-law of index ~ 2.5. Below 4 GeV instead it must be 
significantly suppressed and the total lepton spectrum is dominated by secondary particles. 
The positron spectrum and fraction measured below a few GeV are consistently reproduced 
only within low reacceleration models. We also constrain the scale- height zt of the cosmic-ray 
distribution using three independent (and, in two cases, original) arguments, showing that 
values of < 2 kpc are excluded. This result may have strong implications for particle dark 
matter searches. 
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1 Introduction 

The latest few years have led to remarkable progress in the knowledge of the leptonic compo- 
nent of Galactic Cosmic Rays (CR). This was achieved mainly thanks to a set of successful 
experiments which measured the absolute, combined and relative spectra of electrons and 
positrons in space or in the high atmosphere (see [1] for a recent review). 

Among other relevant results, the PAMELA orbital observatory measured the positron 
fraction between 1 and 100 GeV [2] finding that it unexpectedly increases above 10 GeV; 
Fermi-LAT determined the (e+ + e~) spectrum in the energy range 7 GeV < ii^ < 1 TeV 
[3, 4] to be compatible with a power-law of index —3.08 ± 0.05 and to display significant 
evidence of a spectral hardening above a few hundred GeV; PAMELA recently measured 
the e~ spectrum between 1 and 625 GeV [5] finding that above 30 GeV it is described by a 
power-law of index —3.18 it 0.05 and that at lower energies it is almost in agreement with 
that previously measured by ASM-01 [6]. Furthermore the H.E.S.S. atmospheric Cherenkov 
telescope observed a pronounced steepening of the (e+ + e~) spectrum above 1 TeV [7]. 

All these results can consistently be described in the framework of semi-analytical or 
numerical propagation models accounting for the following spectral components: 

• standard primary electrons (e~) originated by astrophysical sources, most likely Super- 
nova Remnants (SNRs) ; 

• secondary electrons and positrons produced by the scattering of the nuclear component 
of CR with the ISM; 

• a new electron + positron (e^) primary component, of still unknown origin, taking over 
at about a few hundred GeV. 

The extra- component is required in order to provide a consistent description of the rising 
positron fraction found by PAMELA and of the hard (e"*" + e~) spectrum measured by 
Fermi-LAT (see e.g. [8] and Ref.s therein). Its presence has been recently confirmed by the 
measurement of the and e~ separate spectra performed by Fermi-LAT [9, 10]. 



- 1 - 



In spite of these successes several important questions on the origin and the propaga- 
tion of cosmic-ray electrons (CRE) remain open. Due to several degeneracies among the 
relevant parameters, the rigidity dependence of the diffusion coefficient, as well as the role of 
reacceleration and convection in CR propagation are poorly known. 

Another very important source of uncertainty is the large error on the vertical (per- 
pendicular to the Galactic plane) extension of the CR diffusion region (1 < L < 10 kpc) as 
estimated from the ^*^Be/^Be ratio [11]. Interestingly, recent results derived from a global 
analysis of the available data in the framework of numerical propagation yielded a stronger 
constraint L = 5.4 ± 1.4 kpc (corresponding to a 95% posterior range 3.2-8.6 kpc) [12], while 
the same kind of analysis performed with semi-analytic codes found much weaker constraints 
L = 8^7 kpc [13]. We do not discuss here the possible origin of such differences, but just 
remark that the theoretical systematics associated with the determination of this parameter 
may well be larger than the errors estimated by a statistical analysis. Therefore, it is impor- 
tant to use also other, possibly non-local, observables to constrain the scale-height L. These 
issues are relevant not only for conventional CR physics but also for dark matter (DM) indi- 
rect searches since the local flux of DM decay /annihilation products is expected to depend 
significantly on L (see [14] for a recent analysis). 

Furthermore, we know very little about the spatial distribution of CRE sources. In 
particular, one of the most debated questions is whether the source(s) of the extra- 
component is local (e.g. a nearby pulsar) or distributed in a large halo (as expected if its 
source is DM). 

Some of the issues mentioned above can hardly be solved using electrically charged 
messengers alone, because a) solar wind reshapes the charged CR spectra in a poorly known 
way, thereby preventing from probing their local interstellar spectra (LIS) at rigidities below 
a few tens of GeV; h) the diffusion of charged CRs in the random component of the galactic 
magnetic field (GMF) erases most of the angular information about their origin and spatial 
distribution in the Galaxy. 

The 7-ray and the synchrotron diffuse emissions of the Galaxy offer valuable comple- 
mentary probes of the low energy spectrum and of the spatial distribution of CRs in the 
Galaxy. However, the 7-ray diffuse emission receives a significant contribution from proton- 
proton interactions in the interstellar medium. Therefore, we limit ourselves to consider the 
synchrotron emission as it offers a more direct probe of the leptonic CR component. We 
will use the DRAGON numerical CR propagation code to consistently model the electron 
and positron spectra and we will compute the angular distribution and the spectrum of the 
galactic synchrotron emission over a wide frequency range. 

Our approach aims at determining the interstellar spectrum of the CRE below 7 GeV 
under the condition that the Fermi-LAT, PAMELA and H.E.S.S. lepton data are reproduced 
above that energy. We will therefore include the extra-component in our fits. We will 
pay particular attention also to the proton spectrum and to secondary-to-primary ratios by 
requiring that they are all reproduced in our models. This will enable us to constrain relevant 
properties of the source spectrum and of the propagation conditions in the Galaxy. 

The other main goal of this work is to constrain the vertical scale height of the diffusion 
region in the Galaxy. We will show that in addition to the traditional probe ^''Be/^Be, other 
three observables can be exploited and can provide us with valuable insights. Two of them 
yield global (non-local) constraints. The three observables are: 

• the radio spectrum, together with the condition that the GMF strength be compatible 
with Faraday rotation surveys; 
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• the latitude profile of tfie synchrotron emission; 

• the positron spectrum at energies below ~ 5 GeV. 

This paper is structured as follows. In section 2 we will discuss the characteristics of 
galactic magnetic fields. In section 3 we will study the spectra of CRE and of the galactic 
synchrotron emission. In section 4 we will focus instead on the latitude profiles of the galactic 
synchrotron emission and on the problem of the scale height of the turbulent magnetic field 
zt- Finally, in section 5 we will discuss our findings and draw our conclusions. 

2 The Galactic Magnetic Field 

A description of the DRAGON code can be found in [15] and on the DRAGON webpage. ^ 
This code was successfully used to model CR nuclei and lepton propagation in the Galaxy 
and to compute the local interstellar spectra (LIS) of a number of secondary species including 
light nuclei, antiprotons, positrons and 7-rays. 

Differently from other semi-analytical and numerical codes, DRAGON allows for a pos- 
sible spatial dependence of the diffusion coefficient D which we assume to take the general 
form: 



with p = pf3c/{Ze) being the rigidity of the particle of charge Z and momentum p, f{z) de- 
scribes the spatial dependence (in cylindrical coordinates and assuming azimuthal symmetry, 
with 2: = on the Galactic plane and /(O) = 1) of D, and r] provides an effective handling of 
the low energy behavior of D. While one would expect 77 = 1 as the most natural dependence 
of diffusion on the particle speed, several effects may give rise to a different effective behavior. 
Here, where not differently stated, we tune ry, as other parameters, by minimizing the of 
the model against B/C and proton data (see below). 

In previous versions, a simplified single component and z-independent model of the 
GMF was assumed to compute synchrotron energy losses of CR leptons, as that was enough 
to model them with sufficient accuracy. That approximation, however, is not sufficient when 
trying to get a realistic description of the synchrotron emission angular distribution as we 
want to do in this work. Therefore we adopt here a more realistic model of the GMF based 
on Faraday rotation measurements (RMs) of a large number of Galactic and extragalactic 
radio sources (see [16] for a recent review on GMF). 

We consider here two main field components. For the regular component we adopt the 
recent model described in [17] and vary the main parameters in the allowed range (tables 
I and II in that paper). We notice that only the halo component of the regular GMF has 
a role for our analysis, albeit marginal. This has a double torus structure, each of them 
with thickness ~ 0.2 — 0.4 kpc and at distance of 1 — 2 kpc from the Galactic plane (GP). 
Its strength on the torus axis was estimated to be in the range -Bhaio ~ 2 — 12 pG with 
-Bhaio ~ 4 /iG as a best fit. Recently, a rather different model of the regular component of 
the GMF has been proposed [18], which features, besides the spiral arms and disk and halo 
components, also an X-shaped component in the r — z directions, especially relevant close to 
the galactic center. This GMF is explicitly divergence-free. ^ We checked that, as expected, 

^http: //dragon.hepf orge . org/ . 

^The study of the random component associated with this model is in progress [19]. 
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changing from one to the other model does not affect significantly our results. Therefore, we 
will use the model of [17] in the rest of this paper. 

Very little is known on the geometrical structure of the random component of GMF. Its 
energy power spectrum has been determined on the basis of RMs to be [20] 



with a = —0.37 it 0.1 and C = (6.8 it 0.8) x 10"^^ erg cm~^ in the wavenumber range 
0.07 < k < 2 kpc~^, corresponding to a rms field Bq = 6.1 it 0.5 fiG. This measurement 
refers to the Galactic plane region where most of the pulsars are observed. At smaller scales 
the field power spectrum was found to be steeper, compatible with a Kolmogorov one. Its 
normalization matches the extrapolation down to /c ~ 1 kpc of the spectrum measured at 
larger scales. 

Random GMFs are responsible for CR diffusion, hence the spatial and rigidity depen- 
dence of the diffusion coefficient should be related to the turbulent GMF spatial distribution 
and power spectrum, respectively. This is supported by quasi-linear theory and numerical 
simulations of particle propagation in turbulent MFs (see e.g. [21]). 

Concerning the dependence upon rigidity, we remark that it does not need to be the 
same as the one determined from the turbulence spectrum within the quasi-linear theory, 
as noticed in [21]: therefore, we leave the power-law index 6, defined in eq. 2.1, as a free 
parameter. 

We assume instead the same z dependence for i^ran and D~^: 



where zt takes the role of the effective scale-height of the diffusion region (in works by 
other groups a constant diffusion coefficient is assumed, therefore the scale-height zt and the 
position of the vertical boundary of the propagation region L coincide). In order to ensure 
that boundary effects are not relevant, we solve the transport equation in DRAGON by setting 
always the vertical boundary L = Szt- 

3 Cosmic-ray electron, positron and synchrotron spectra 

In this section we use DRAGON to model the CR electron and positron spectra and to 
consistently compute the synchrotron emission of the Galaxy integrated along the line of 
sight. 

In order to assess the uncertainties related to CR propagation, we consider four repre- 
sentative propagation regimes: PD (plain diffusion), KRA (Kraichnan), KOL (Kolmogorov) 
and CON (convective) . For each of them we vary the scale-height of the diffusive halo in the 
range = 1 16 kpc. The main parameters of these setups for zt < 8 kpc are summarized 
in table 1 . All these models are chosen in such a way as to minimize the combined against 
the B/C and the proton observed spectra, as shown in figure 1. Their reduced is always 
smaller than unity. This approach is similar to the one adopted in [14], to which we refer the 
reader for more details including the used data set. 

Once the model parameters are fixed according to this procedure, we have still the 
freedom to choose the spectral index and normalization of the injection spectrum of the 
primary electrons and of the extra-component. We fix these parameters by fitting the e~ +e'^ 




(2.2) 



Diz) OC -Bran(z) OC exp{-z/zj} 



(2.3) 
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Figure 1. B/C and proton spectra computed for the reference propagation setups PD4 (black lines), 
KRA4 (blue) K0L4 (red) and C0N4 (purple) defined in table 1 are compared with a selection of 
experimental data. The modulation potential used for the B/C is $ = 0.5 GV. 



spectrum and the positron fraction measured by Fermi-LAT. In contrast, the spectrum and 
normalization of the secondary electron and positron spectra are determined once the proton 
and helium spectra (and the propagation parameters) are fixed. 

Solar modulation is one of the major issues in modeling CR spectra below ~ 20 GeV. 
The most common treatment assumes the force-field charge independent approximation [22] . 
In this approach the effects of propagation in the solar system on the observed spectra 
of all electrically charged CR species are described in terms of the single, time-dependent 
parameter $ (the so called modulation potential), to be fixed against CR data. However, 
charge dependent drifts in the complex magnetic structure of the heliosphere may lead to 
significant deviations from this simple scenario (see e.g. [23, 24] for recent discussions of this 
problem) . 

In order for our analysis to be as independent of solar effects as possible, we tune 
the CRE source parameters by comparing the force field modulated e~ + e"*" spectra with 
Fermi-LAT data, which have been measured above 7 GeV. According to [25], and taking 
into account that data were taken during a quiet Sun phase, drifts effects are expected to 
be small. Below that energy we will model the e~ + e"*" LIS on the basis of the observed 
synchrotron spectrum of the Galaxy which is unaffected by propagation in the heliosphere. 
At this point the low energy CR electron and positron spectra at Earth will be predicted 
(assuming a model for solar effects) and low energy data can be used to further constrain 
the model. 

A comment is in order at this point: We could have chosen to consider the electron 
spectrum and the positron fraction measured by PAMELA in order to fix the relevant CRE 
parameters, and then compare with the synchrotron spectrum and the Fermi CRE spectra. 
We did not follow this route because: a) PAMELA data on the CRE spectrum do not extend 
to high enough energy to probe accurately the extra-component; b) PAMELA data extend 
instead down to low energies, where solar modulation effects can be significant. We remark 
that using a combination of the spectra observed by both experiments is not feasible either, 
because of the presence of a systematic offset between PAMELA and Fermi data, which 
would make extremely difficult a proper comparison between them. 

After having computed with DRAGON the CR electron and positron spatial distributions 
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iVlouel 


Do 

(1028 cmVs) 


(km/s) 


dVc/dz 
(km/s/kpc) 


6 


rj 


7 


% 
(GV) 


PD8 


3.64 


0* 


0* 


0.56 


-0.42 


2.3 


0.57 


PD4 


2.24 


0* 


0* 


0.57 


-0.40 


2.29 


0.56 


PD2 


1.11 


0* 


0* 


0.58 


-0.45 


2.28 


0.57 


PDl 


0.515 


0* 


0* 


0.58 


-0.40 


2.29 


0.58 


KRA8 


4.28 


15.7 


0* 


0.5* 


-0.39 


2.35 


0.68 


KRA4 


2.60 


13.7 


0* 


0.5* 


-0.39 


2.36 


0.65 


KRA2 


1.35 


15.6 


0* 


0.5* 


-0.41 


2.35 


0.67 


KRAI 


0.627 


15.3 


0* 


0.5* 


-0.40 


2.33 


0.69 


K0L8 


7.17 


38.9 


0* 


0.33* 


1* 


2.00/2.40 


0.52 


K0L4 


3.98 


33.4 


0* 


0.33* 


1* 


1.93/2.47 


0.51 


K0L2 


2.07 


31.9 


0* 


0.33* 


1* 


1.91/2.47 


0.54 


KOLl 


1.06 


35.9 


0* 


0.33* 


1* 


1.81/2.40 


0.47 


CONS 


0.952 


36.1 


50* 


0.58 


1* 


1.92/2.48 


0.52 


C0N4 


0.923 


36.0 


50* 


0.53 


1* 


1.94/2.48 


0.52 


C0N2 


0.794 


34.5 


50* 


0.48 


1* 


1.95/2.46 


0.52 


CONl 


0.573 


31.4 


50* 


0.41 


1* 


1.95/2.46 


0.52 



Table 1. We report here the main parameters of the reference propagation setups used in this work. 
The number appearing in the name of the model indicates the used value of zt- The convective velocity 
Vc is assumed to vanish on the Galactic plane (z = 0). The KOL and CON models have a break in 
the nuclei source spectra 7 at a rigidity of 11 GV. Parameters marked with a * are fixed a priori to 
characterize the setup. More details on these models are reported on the DRAGON webpage. 

and energy spectra in the Galaxy, we use the standard formalism to compute the synchrotron 
emissivity at each point. For a monochromatic and isotropic distribution of electrons, this 
reads [26, 27] 



e(l^,7)rcg 
e(j^, 7)rand 



V3- 



Cy 



mc 
2 



2 -^pcrp-^(^) 



(3.1) 
(3.2) 



In eq. 3.1, e(z^, 7)i.og refers to the unpolarized emissivity at the frequency v due to the regular 
magnetic field component perpendicular to the line of sight -Bperpj from electrons with energy 
E = mc^^. We define x = v/v^c^, where i/c^ = 3e/(47rmc)i?perp7^ is the critical synchrotron 
frequency, while F{x) = x K^i^{x')dx' . The functions Ki are modified Bessel functions 
or order I. In the expression of the isotropic emissivity due to the random component we 
have C = 2\/3e^/(mc2) Brand, where y = iz/i/^"""^ and z^^'^'''^ = 3e/(27rmc)-Bi.and7^- For each 
line of sight we then integrate the emissivity to obtain the intensity 

Sometimes, we will refer to the brightness temperature, which is defined as 



T 



I{vy 



B 



(3.3) 



where ks is the Boltzmann's constant. 

In the rest of this section we will keep the diffusive halo scale-height fixed at the reference 
value Zt = 4: kpc. The effect of varying this parameter will be discussed in section 4. 
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3.1 CRE models based on Fermi-LAT and synchrotron spectrum data 

As we mentioned in the Introduction (see also [4, 8, 10]), an acceptable fit can be found only if, 
besides a standard electron e~ component with a power-law source spectrum Jq oc E~"'°^^ \ 
we introduce an electron and positron charge symmetric extra-component with the form 
JS (E) oc E'^^^^exp{-E/Ecnt) with ^cut ^ 1 TeV. We assume that the spatial distribu- 
tion of both the standard and extra-component sources traces that of SNRs as determined 
from pulsar catalogues. In section 5 we will briefly discuss how our results change assuming 
the extra-component to be local hence not contributing to the large scale Galactic syn- 
chrotron emission. We also account for secondary electrons and positrons produced by CR 
nuclei scattering onto the ISM (only p — p, p — ^He and ^He — ^He are relevant) adopting 
the parameterization of the production cross sections provided in [28]. The cross sections 
provided in [29] do not differ from the ones in [28] in the energy range relevant to our work 
[30]. 

We then compute the synchrotron spectrum of the Galaxy for the CRE models obtained 
as described above. For each of these models we tune consistently, through an iterative pro- 
cedure, the normalization of the random component of the GMF in order to fit the observed 
synchrotron spectrum at 408 MHz where radio data are most complete and a possible con- 
tamination due to free-free emission is negligible [31]. For the regular component of the GMF 
we adopt model [17] (see section 2) which was built on Faraday rotation measurements. In 
figure 2 we show for the KRA4 that the contribution of that field component to the syn- 
chrotron spectrum is subdominant. The same holds for the other setups. As a consequence, 
the effect of changing the regular GMF model among those discussed in section 2 is negligible 
as we checked also for different values of zt- 

We integrate the Galactic emission along the line of sight and use HEALPix ^ to properly 
average the resulting flux over the sky region 40° < I < 340°, 10° < b < 45° where / and b are 
Galactic longitude and latitude respectively. This is the region where the contamination from 
point-like and local sources is expected to be the smallest and it is the same considered in [32] . 
In this region we compare the simulated spectra with the ones measured by a wide set of radio 
surveys at 22, 45, 408, 1420 , 2326 MHz as well as WMAP foregrounds at 23, 33, 41, 61 and 
94 GHz as consistently catalogued in [31] (to which refer the reader for the survey references) 
where they were pixelized and transformed in Galactic coordinates. Above 1 GHz the free-free 
emission was estimated to give a non-vanishing contribution which however is subdominant in 
the considered frequency and latitude intervals; the 94 GHz channel is likely to be dominated 
by spinning dust-emission. Absorption due to free-free scattering is estimated to be negligible 
down to 10 MHz [33]. The error bars on the radio fluxes reported in our figures correspond 
to the flux semi-dispersion in the considered sky region. We remove the Galactic plane, as 
well as discrete and local diffuse radio sources (e.g. the Galactic Spur) using the extended 
temperature analysis mask used by the WMAP collaboration. ^ The CRE models obtained 
according to this procedure are shown in table 2. 

Figure 2 (solid lines) shows a set of representative models which fit Fermi-LAT data and 
provide simultaneously a good description of the measured synchrotron spectrum. Noticeably, 
for all setups with = 4 kpc the required strength of i?ran(0) is well within the interval 
allowed by RM (see section 4 for a wider discussion on this issue). 

We notice, however, that the experimental data are reproduced only if a sharp spectral 

■''http: //healpix . jpl .nasa.gov. 

^http: //lambda. gsfc .nasa.gov/product/map/dr4/masks_get . cfm. 
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Model 


70 (e ) 


7(e±) 


^cut (GeV) 


$ (GV) 


Xfermi 


PD4 


0.9/2.53 


1.38 


986 


0.3 


0.128 


KRA4 


0.6/2.53 


1.50 


967 


0.3 


0.20 


K0L4 


0.5/2.54 


1.56 


985 


0.1 


0.177 


C0N4 


0.5/2.45 


1.45 


999 


0.1 


0.274 



Table 2. We report here the source parameters which yield a consistent fit of the CRE and 
synchrotron spectra. The reported values of low energy spectral indexes are the maximum allowed 
to reproduce radio data at 2ct. The values of the force field modulation potential $ have been fixed 
to match the e~ AMS-01 data. The reduced is computed against the e+ + spectrum and the 
positron fraction measured by Fermi-LAT. 




Figure 2. The average synchrotron spectra in the region 40° < I < 340°, 10° < b < 45° computed 
for the reference propagation setups PD4 (black lines) , KRA4 (blue) K0L4 (red) and C0N4 (purple) 
defined in table 1 are compared with experimental data derived from [31]. For each setup we show the 
spectra obtained with (solid lines) and without (dashed line) the spectral break in the e~ source spectra 
given in table 2. Dotted lines represent the corresponding contribution of secondary e~ and e"*". The 
contribution to the synchrotron fiux of the regular GMF with Bhaio = 4 /iG (see [17]), computed for 
the KRA4 setup, is shown as the triple dot-dashed line. The random component is tuned to reproduce 
the spectrum normalization at 408 MHz. The required normalization is -Bran(O) = 7.6 /xG. The critical 
synchrotron frequencies, calculated for this value of i3ran(0), are reported for a few reference values 
of the electron energy. 



hardening, a break, is imposed to the e source spectrum below ~ 4 GeV. Its amplitude 
must be A7(e") ~ 1.6(1.9) for the PD (KRA) setups and A'y{e") ~ 2.0 for the KOL and 
CON (see table 2). This is similar to what was found in [32] with Aj{e~) = 0.9, and in 
[34] where A'j{e~) = 1.57. We find some significant differences with respect to these papers 
though. 

In contrast to [34] we find the required break to be significantly dependent on the 
propagation setup (see table 2). However, in [34] the break of the propagated spectrum was 
quoted, rather than that of the injected spectrum. It is quite natural that the way in which 
the break translates into the other depends on the propagation model. Moreover, the break 
quoted in [34] is that of the total (primary+secondary) propagated spectrum. However, 
the contribution of secondary electrons and positrons depends on the primary proton and 
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helium spectra in addition to the propagation setup (see the dotted lines in figures 2 and 3), 
and cannot be arbitrarily altered. Similarly to what we noticed in [35] (see also [36] for an 
analogous remark), the spectra of secondary particles display a bump at £^ ~ 1 GeV. The 
bump is particularly prominent for models with strong reacceleration, as the KRA and the 
CON setups. Due to the presence of this component, the observed synchrotron and CRE 
spectra can be reproduced only by imposing a strong break or a cutoff in the primary e~ 
component. 

Differently from [32] we find a larger flux of secondary particles in the GeV range and 
that a more pronounced break is required. We notice that the values of the low-energy 
spectral indexes reported in table 2 and adopted in figure 2 (solid lines), are maximal values, 
i.e. they correspond to the minimal breaks, allowed by the radio data at 2a. Indeed we find 
that assuming lower values of 70 (e~), or replacing the break with an exponential infrared 
cutoff of the primary e~ component at -Eiow ~ 2 GeV, do not change our results significantly. 
Therefore, below ~ 100 MHz the synchrotron spectrum of the Galaxy becomes gradually 
dominated by secondary e~ and e"*". 

3.2 Back to the e+ and e~ observed spectra 

In the last part of this section we turn our attention back to the electron and positron spectra 
measured in the solar system. Our aim is to check to which extent the models we built to 
describe the synchrotron emission are consistent with CRE data even below 7 GeV. This 
results will provide a valuable consistency check for the validity of our approach and of the 
force field modulation treatment that we use in this paper. 

We only consider data taken only during periods of low solar activity since under those 
conditions force field modulation is expected to provide the most reliable results. Our ref- 
erence data sets below 7 GeV are those taken by AMS-01 [6] in June 1998 and PAMELA 
between July 2006 and January 2010, both during low activity periods. We notice that the 
latter period is different from that during which PAMELA measured the proton spectrum 
which we used to tune our propagation setups (2006-08). This explains why the modulation 
potentials reported in tables 1 and 2 are different. 

For each propagation setup, Fermi-LAT data can be fitted for several combinations 
of the e~ source spectral index and the modulation potential Low energy electron and 
positron data, however, break this degeneracy. The reason is that, as follows from our 
previous considerations, in the GeV and sub-GeV regions the e~ and e"*" fluxes receive large, 
or dominant, contributions from secondary particles which are unaffected by the primary 
source spectrum and, in addition, the total spectrum is fixed by the synchrotron spectrum. 
The modulation potential is then fixed against AMS-01 data. In conclusion, for each setup 
the models in table 2 are unambiguously fixed by the requirement to match the observed CRE 
and synchrotron spectra. Figure 3 show that all models nicely reproduce the e~ spectrum 
measured by AMS-01 down to ~ 1 GeV. PAMELA e~ data are also reproduced if we allow 
for a normalization offset compatible with the systematic errors of PAMELA and Fermi-LAT 
(relative to which our models are normalized). 

Our findings for the positron flux are even more interesting. The spectrum below 
10 GeV strongly depends on the propagation setup. The KOL and CON configurations are 
clearly incompatible with AMS-01 data because of the presence of a pronounced bump at 
E ^ \ GeV, arising due to the strong reacceleration (see also [35]). The PD4 and KRA4 
models, on the contrary, match the data more closely. 
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Figure 3. The + e~ (panel a), e~ (panel b) and 6+ (panel c) spectra multiplied by as well 
as the positron fraction (panel d) computed for the reference propagation setup PD4 (black lines), 
KRA4 (blue), K0L4 (red) and C0N4 (purple) defined by the parameters in tables 1 and 2 are shown 
together with a selection of experimental data sets. Continuos and dashed lines represent modulated 
and interstellar (LIS) spectra respectively. The dotted lines correspond to the secondary contributions 
to the LIS spectra. 

For what concerns the positron fraction, there is a clear (and well-known) discrepancy 
between AMS-01 and PAMELA data below a few GeV. While this may be partially due 
to a normalization offset between the two experiments, the different slopes at low energies 
can only be explained by solar modulation and the drift of the particles in the heliosphere 
[37]. Indeed, these experiments were operating during periods with opposite polarity of the 
heliospheric magnetic field. It is a quite generic feature of drifts in the solar system to either 
keep particles diffusing tight to the solar current sheet (when the charge sign and the polarity 
are opposite) or to push them outside the current sheet (when the charge and the polarity 
have the same sign) [23, 24]. This will therefore increase or decrease the propagation time in 
the heliosphere, hence enhancing or lowering the energy losses, respectively. 

When dealing with ratios like the positron fraction, the situation is even more compli- 
cate, because the effects described above mix (an account of these arguments related to the 
very similar case of the antiproton/proton ratio can be found in [38] and references therein). 



- 10 - 



Interestingly, it was found that in phases of low solar activity and positive polarity, as that 
when AMS-01 was operating, the positron fraction computed with the force field approxima- 
tion and a small modulation potential agrees quite well with those models [37]. Our models 
are based on high energy CRE data and synchrotron spectra, hence they are not affected 
by drift effects in the solar system. Therefore, the agreement with AMS-01 data down to 
~ 1 GeV provides a valuable check of the consistency of our approach. 

We notice that the KRA model in [39], though similar to its homologous in this paper, 
reproduced the positron fraction measured by PAMELA down to 1 GeV rather than that of 
AMS-01. That model, however, adopted a weaker break of the e~ source spectrum as the 
synchrotron emission of the Galaxy was not fitted there. This shows the relevance of the 
latter channel to understand low energy CRE spectra. 

We conclude that low reacceleration models with a strong flattening, or a IR cutoff, in 
the e~ injection spectrum below a few GeV are consistent with the full set of radio and CR 
direct measurements. 

4 Halo height and the latitude profile of the synchrotron emission 

4.1 The effects of varying Zt on the synchrotron spectrum and angular distri- 
bution 

We study here how our previous results depend on the scale height zt of the GMF random 
component and place constraints on this quantity on the basis of available experimental data. 

Our first argument exploits the Zf-dependence of the synchrotron emission flux at 
medium and high Galactic latitudes. The flux observed in a given direction is proportional 
to the line of sight integral J B'^{l)ne{l)dl, where Ue is the e~ + e"*" CR density and the inte- 
gral extends to the boundaries of the computational domain. In figure 4 we show the ne{z) 
profile computed with DRAGON for several propagation setups and different values of zt- At 
~ 1 GeV energy losses have a minor effect on the CRE large scale distribution. In fact, at 
that energy they are dominated by IC scattering (with the exception of the GP region where 
bremsstrahlung dominates). The IC diffusive loss length is of the order of a few kpc on the 
GP and increases rapidly away from the disk. We showed in the previous section, that i?ran 
is the dominant GMF component contributing to the synchrotron emissivity. Consequently, 
for a given propagation setup, the synchrotron fiux depends only on the random field nor- 
malization i?ran(0) and on zt- Measurements of the former quantity which are independent 
on the synchrotron emission offer, therefore, a valuable probe of the latter. 

We consider then the synchrotron spectrum. For each value of zt we tune -Bran(O) so 
that the computed spectrum reproduces the data at 408 MHz. Interestingly, we find a tight 
relation i?^g^^(0) oc z^^ (see figure 5). This relation corresponds to the physical requirement 
that the total energy of the GMF be independent of the choice of zt- In other words, the 
synchrotron spectrum considered here is sensitive to, and fixes, the total magnetic energy of 
the galaxy (see [40] for a similar discussion). 

In figure 5 we report the 3 and 5 a regions allowed by RMs of Galactic pulsars [20] 
(see section 2). Values of zt smaller than 4(3) kpc are excluded at 3 (5) a. Thick CR halos, 
however, cannot be excluded on the basis of this argument. The plot in figure 5 is computed 
for the KRA setup. We find no significant differences between this setup and the PD so that 
the same constraint applies to different propagation conditions. 

A complementary probe of the vertical distribution of CRE is offered by the latitude 
profile of the synchrotron emission. In figure 6 we compare the observed latitudinal profile 



- 11 - 




Figure 4. Top panels: The vertical profile of the e~ density at i? = 4 kpc is plotted for E = 1 
and 10 GeV and for different propagation setups. The C0NST4 (green lines) setup is similar to the 
KRA4 (and it also fit B/C and proton data), but the choice of a uniform diffusion coefficient and a 
random GMF up to ^max = Zt- Bottom panel: the same quantity computed for the KRA setup with 
zt = 1,2,4,8, 16 kpc at E = 1 GeV . Fluxes are normalized to their values on the galactic plane at 
z = 0. 

of the synchrotron emission at 408 MHz to that calculated for the KRA setups for different 
values of zt. As done before, we look at the region 40° < / < 340°, 10° < 6 < 45°. Error bars 
correspond to the semi-dispersion of the distribution of the observed flux in those regions. 
We also consider the subregion 40° < I < 100°, 10° < b < 45°, where the flux variance 
(which is likely dominated by non-subtracted local structures) is smaller, which yields more 
stringent constraints. For each zt we tune the value of -Bran(O) so that the average spectrum 
in these regions is reproduced. We see that low values of zt are disfavored. This agrees with 
the results of the spectral analysis. A analysis shows that zt < 2 kpc are excluded at the 
3a level. The relevant values are reported in table 3. 

One of the most serious theoretical uncertainties which may affect the results of this 
subsection arises from the poorly known form of the vertical profile of the random GMF, hence 
of the diffusion coefficient. In order to quantify this uncertainty we repeated our analysis using 
a Gaussian profile (rather than the exponential used above) and, as an extreme possibility, 
a step function profile f{z) = 9{\z\ — zt) (which is the profile adopted in [32, 34]) for both 
quantities. As a consequence, smaller field strengths are required on the GP to reproduce the 
synchrotron spectrum normalization and the curve in figure 5 translates downwards. This 
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Figure 5. The normalization of tlic random GMF is plotted against its vertieal scale height (whieh 
we assume to be the same as that of the diffusion coefficient). The 3(5) a regions allowed by RM 
data arc represented in gray (light-gray). Red dots are our results obtained under the condition that 
KRA models reproduce the observed synchrotron spectrum. The black line is a B^an oc f/zt fit of 
those points. The fitting function is (-Bian/1 mG)'^ = 148.06 (1 kpc/zt) + 19.12. The fit computed for 
the other setups considered in this work would superimpose to that line. 
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Figure 6. The latitude profiles of the synchrotron emission at 408 MHz in the region 40° < I < 100° 
computed for the KRA propagation setup and zt = 1,2,4,8,16 kpc are compared with radio data. 
The grey shadowed region is not considered when placing the constraint. 

would turn into a smaller lower limit on zt, from 3 to 2(1) kpc at 5 o" for the Gaussian 
(step-like) profile. We also find that zt > 9(10) kpc are excluded at 3cr (but still allowed 
at 5a). We notice, however, that a step-like profile is physically hard to justify, as the gas 
density, which hosts the currents that should sustain the magnetic fields, decreases rapidly 
outside the galactic plane. 
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zt (kpc) 




1 


2.7 


2 


2.5 


4 


1.6 


8 


0.9 


16 


0.4 



Table 3. Reduced values (13 degrees of freedom) of the comparison between synchrotron profiles 
and observational data in the region 40° < I < 100°, 10° < 6 < 45°. 
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Figure 7. The latitude profiles of the synchrotron emission at 408 MHz in the regions 40° < I < 100° 
computed for the KRA4 propagation setup with an exponential profile of the random component 
of the GMF (blue line) is compared with that obtained with similar models but a step-like profile 
(C0NST4) (green line) and a Gaussian profile (red line). The grey shadowed region is the same as in 
figure 6. 

We also consider the effect of changing f{z) on the latitude profile of the synchrotron 
emission. In Fig. 4 we already compared the CRE profiles computed for the KRA4 with the 
C0NST4, which adopts a step-like f{z), finding small differences. In figure 7 we see as the 
latitude profile of the synchrotron emission obtained for those setups differ by 20% at most 
(the Gaussian profile which is also shown in that figure is even closer to the exponential one 
used for the KRA4 setups). Therefore, we conclude, that observable offers a reliable probe 
of the CR diffusive halo scale height. 

4.2 The effects of varying Zt on the electron and positron spectra 

In section 3.2 we showed that the low energy CR positron spectrum can be conveniently 
used to probe some of the CR propagation properties, in particular the effectiveness of 
reacceleration. Here we demonstrate that the same can be done also for zt- 

We remind the reader that the B/C, or other secondary-to-primary ratios, can constrain 
only the ratio Do/zt, where Dq is the normalization of the diffusion coefficient (see eq. 2.1). 
The secondary flux however depends significantly on the scale height zt, thus breaking 
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the degeneracy of this parameter with Dq. This can be promptly understood from the simple 
general equation (see e.g. [41]) 

iV(e) = ^ (4.1) 

where r is the minimum between the escape and loss times, Q the production rate per unit 
volume, Vs and V the source and CR filling volumes respectively. For secondary electrons 
and positrons Q oc CTppHgasNp where Np is the CR proton density in the GP. We showed 
that at the most relevant energy ~ 1 GeV the CRE fill the entire volume occupied by 
the random GMF, so that V oc zt- This also implies that energy losses are subdominat and 
that r is determined by the propagation setup through the D^/zt ratio (thus Tesc)- As a 
consequence, for a given setup the secondary electron and positron densities decrease with 
increasing zt- Since this is a dilution effect, it is expected to be energy independent which 
is what we find from our numerical computations (see figure 8). This effect is most clearly 
detectable in the e"*" channel since below 10 GeV the contribution of primary positrons from 
the extra-component is negligible. Figure 8 shows that the e+ spectrum measured by AMS- 
01 disfavors values of zt smaller than 2 kpc. This strengthens the conclusions reached in the 
first part of this section on the basis of radio observations. In figure 8 we only show spectra 
computed with the KRA models since the other setups yield similar results. Remarkably, 
even large values of zt cannot reconcile strong reacceleration models with the positron data, 
because the dilution effect cannot compensate for their typical bump at ~ 1 GeV. 



5 Discussion and conclusions 



One of the main implications of this work is that below a few GeV the injected electron 
spectrum must be drastically suppressed with respect to a simple power-law extrapolation of 
the one which fits CRE Fermi-LAT data. This suppression can be obtained by introducing 
either a break or an exponential IR cutoff in the e~ source spectrum below a few GeV. This 
last possibility is particularly intriguing, as it might hint at a suppression of the leakage of 
accelerated electrons from the sources at GeV and sub-GeV energies. Furthermore, similarly 
to [32], we find a rather high source spectral index 70 (e~) ~ 2.5 above a few GeV. This is 
larger than the value inferred from SNR radio catalogues derived in [42] which is 7snr = 
2.0 lb 0.3, albeit the two are compatible within 2a. As discussed in [32], this may imply that 
the electron population probed by the SNR radio observations is not representative of that 
injected in the ISM. However, we point out that this has to be considered as an effective 
index, since the code does not take into account the presence of structures in the source 
distribution. Accounting for the spiral arm distribution of SNRs may result in a different 
requirement for the injection index: indeed, being the Sun in the so-called local spur situated 
in a interarm region, the average distance from SNRs is larger than in the smooth case: as 
a consequence, a harder injection spectrum may be required to compensate for the larger 
energy losses and reproduce the observed spectrum. 

We find a break amplitude similar to that in [34]. In that paper, however, the break 
was imposed on the total (primary + secondary) propagated spectrum. It was claimed 
that leaving the secondary spectra unaffected would lead to overproduction of synchrotron 
radiation. Our figure 2, however, shows that this is not the case. Below ~ 100 MHz the 
contribution of secondary electrons and positrons becomes gradually sufficient to account for 
the observed radio spectrum and no break is required in the spectra of these particles. For 
opposite reasons, we disagree with the results of [32], as in that paper it was found that 
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Figure 8. The e"*" + e~ (panel a), e~ (panel b) and e"*" (panel c) spectra multiplied by as well as 
the positron fraction (panel d) spectra computed for the KRA setup and zt = 1, 2, 4, 8 kpc (from top 
to bottom respectively) are shown against a selection of experimental data. The line style notation is 
the same as in figure 3. All models are modulated with $ = 0.3 GV. 

secondary leptons can contribute at most one-third of the low energy synchrotron spectrum. 
This difference can be explained by the break in the proton and Hehum injection spectra 
in their models [43], which suppresses secondary production at low energy. We remark that 
we adopt a break only in KOL and CON models, while in the other, experimentally favorite 
models the primary nuclei spectra are a single power-law in momentum. 

A relevant consequence of our result is that an interpretation of the low energy (low 
frequency) behavior of the CRE (synchrotron) spectrum in terms of a break in the rigidity 
dependence of the diffusion coefficient (which should affect all components) is disfavored. 
A strong break, as proposed in [44], is hardly compatible with both CR nuclei and lepton 
spectra (a multichannel analysis is not provided in that paper). 

One relevant aspect of our approach is that it is based on a multichannel analysis of 
nuclear, electron and positron spectra. Positrons are especially important because their low 
energy spectrum is dominated by secondary particles. We find that once the low energy 
e~ source spectrum is tuned to reproduce the observed synchrotron spectrum, only models 
featuring low reacceleration can reproduce the observed e"*" spectrum and fraction. This 
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conclusion is in agreement with [32] and it strengthens the results of our previous analyses 
based on antiprotons [39] and e"*" + e~ spectra [35]. 

Our modeling of the synchrotron emission of the Galaxy accounts, for the first time 
in this framework, for the presence of the extra-component. This is required not only 
to consistently model PAMELA and Fermi-LAT high energy data (E ^> 10 GeV) but also 
to correctly estimate the e~ source spectrum from CR and radio data. In this work we 
treated the extra-component source term as a continuous distribution and assumed it to 
spatially coincide with that of the standard e~ component. This may be not realistic, because 
the extra-component could well be powered by one or a few nearby sources [8]. We find 
however that adopting the latter hypothesis will not affect significantly any of our results. 
Indeed, we have checked that subtracting the contribution of the extra-component from the 
-|- e~ Galactic spectrum will not change the synchrotron spectrum at a detectable level. 
However, by combining high precision CRE data from AMS-02 [45] and CALET [46] as well 
as forthcoming radio observations by PLANCK^, LOFAR^ and SKA'' can become possible 
to detect the contribution of the extra-component to the synchrotron spectrum. 

With our method we exploit the galactic diffuse synchrotron emission as a handle to 
measure the low energy LIS spectrum of CR electrons and positrons, in a similar way as 
the diffuse 7-ray emission is being exploited to obtain the LIS proton spectrum [47-49]. 
This is a valuable information for studies of solar modulation. While the simple force-field 
approach works quite well in reproducing AMS-01 data, PAMELA data require a more 
detailed investigation of charge-sign dependent effects in the solar system. 

We investigate the vertical extension of the CRE distribution by means of several com- 
plementary methods. Remarkably, these methods are more powerful than the ^''Be/^Be in 
placing lower limits on zt and also more robust because: 

• They come from non-local observables and they are not sensitive to the effects of solar 
modulation or local physics (e.g. the possible presence of a local bubble [50]); 

• the constraints based on ^'^Be/^Be strongly depend on the propagation model. 

Given the importance of these points, we show for reference in figure 9 the vertical scale 
dependence of the ^'^Be/^Be in our KRA and CON models, together with available experi- 
mental data from ACE and ISOMAX. The most discriminatory points are from ACE where 
they are susceptible to solar modulation. The dependence on the propagation setup is clearly 
seen by comparing the KRA and the CON model. Qualitatively, the difference between the 
two cases is that in convective models the convective outflows along the normal to the disc 
plane introduce an effective vertical scale height [51, 52]. This explains why in the CON case 
the ^'^Be/^Be is less sensitive to zt than in the KRA case. As a consequence, the constraining 
power of ^''Be/^Be is reduced for convective models. The results we obtained from our new 
analysis are not plagued by such uncertainties. 

Having lower limits on the halo scale height is of great importance for indirect DM 
searches. Because the sources of DM originated particles are distributed also in the halo (in 
contrast with astrophysical sources that are more localized in the disk), the predicted fluxes 
depend strongly on the extension of the halo. In particular, constraints on the annihilation 
cross section are considerably weakened if the halo size is small. According to [14], imposing 

^http : / /planck .esa.it 
®http: //www . lof ar . org 
''http : / /www . skatelescope . org 
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Figure 9. ^'^Bc/^Be modulated (dotted) and non-modulated (solid) spectra for the KRA (left panel) 
and CON (right panel) configurations. From bottom to top curves correspond to Zt — 1, 2, 4, 8 kpc. 



the limit zt > 0.5 kpc would improve the constraints on the annihilation cross section by a 
factor of two to ten. A full account of the consequences of our bound zt > 2 kpc is beyond 
the scope of this paper and will be left for a forthcoming work. 

Our lower limit on zt may have also important implications for the interpretation of the 
CR anisotropy. In fact, it was shown that increasing the diffusive halos scale height reduces 
the possible role of source stochasticity to explain the large scale anisotropy fluctuations 
observed above 10 TeV [53]. 

For the first time, we have placed a constraint on the CR diffusive halo scale height 
based on the comparison of the computed synchrotron emission intensity with observations. 
It is based on the CRE spectra measured by Fermi-LAT and GMF constraints determined 
with Faraday RMs. Our constraint depends weakly on the poorly known vertical profile 
of the random component of the GMF being zt < 3(2) kpc for a exponential (Gaussian) 
profile. Statistical and systematic errors on i?ran are expected to be significantly reduced in 
the near future by LOFAR and SKA. Similarly to what was done in [34] we also compare the 
computed latitude profile with radio maps of selected regions, obtaining consistent results 
and excluding < 2 kpc at 3cr confidence level. We remark that constraints derived in this 
way are stronger than the ones obtained from 7-ray data [? ]. 

Finally, we used the low energy positron spectrum and positron fraction dependence on 
Zt to constraint this quantity. Also this method seriously disfavors low values of zt- AMS-02 
data and recent progresses in solar modulation modeling should soon considerably improve 
the usefulness of this approach. 
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